A new basis set for the description of electrons in superconductors 
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In the usual description of electrons in metals and superconductors, the single electron states are 
assumed to satisfy Bloch's theorem. This is because the electron-ion interaction is privileged over 
the electron-electron interaction. However the theory of hole superconductivity proposes that in 
the transition to superconductivity carriers 'undress' from the electron-ion interaction. I propose 
here a new basis set to describe electrons in the superconducting state that does not satisfy Bloch's 
theorem but instead is designed to optimize the electron-electron interaction. The new basis set 
favors a superconducting state where a spin current exists and where states near the bottom of the 
band become partially occupied, as predicted by the theory of hole superconductivity. 
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I. INTRODUCTION 

It is generally assumed that electron states in a peri- 
odic solid obey Bloch's thcorcni^l]. That is, a wavevector 
k characterizes the single electron orbital 'i/'nk(r) which 
obeys the relation 



e*-^^„k(r) 



(1) 



where the index n (band index) differentiates different 
states with the same Bloch wavevector k. In Eq. (1), R 
is any lattice translation vector that leaves the electron- 
ion potential invariant: 



f/e_.(r + R) = (r). 



(2) 



Electron-electron interactions are assumed to convert the 
single-particle electron-ion potential to a new 'effective' 
potential Uef / that includes electron-electron interaction 
direct, exchange and correlation effects (e.g. the Kohn- 
Sham potential in density- functional theory) Q. The ef- 
fective potential still obeys 



C^e//(r + R)-C/e//(r) 



(3) 



with the same lattice vectors R, and the corresponding 
single particle orbitals (e.g. Kohn-Sham orbitals) still 
obey Bloch's theorem Eq. (1). 

Similarly, in the generally accepted conventional the- 
ory of superconductivity (BCS theory) Q the many-body 
wavefunction is written as 



I* >Bcs^ Yliuk + "k4^c^ki)|0 > 



(4) 



where the index k is the Bloch wavevector of the state 
where the single electron of spin a is created by the opera- 
tor cj^^. The normal state of the system can be expressed 
in the form Eq. (4) with itk^k = and Wk = or Uk = 
depending on whether ek < or ek > ej? with ep the 
Fermi energy. The superconducting state instead is char- 
acterized by UkWk 7^ for some values of k. For the case 



of nearly free electron metals, the Bloch states are often 
simply taken to be free electron plane wave states. 

However, it is not necessarily true that electron states 
in a periodic electron-ion potential have to obey Bloch's 
theorem. Consider a case where the electron-electron in- 
teraction completely overwhelms the electron-ion interac- 
tion. The electronic spatial arrangement that minimizes 
the energy may be an electronic lattice (Wigner crys- 
tal, incommensurate charge density wave) with a lattice 
structure that has no relation with the underlying ionic 
structure. For example, the electronic structure may be 
hexagonal and the underlying ionic lattice cubic. In such 
situation the single particle states describing the elec- 
trons will not obey Bloch's theorem Eq. (1). 

We have recently argued that experimental evidence 
clearly indicates that in the transition to superconduc- 
tivity the carriers undergo a wavelength expansion and 
end up completely 'undressed' from the electron-ion in- 
teraction, so that they no longer 'see' the discrete nature 
of the electron- ion potential^, It is natural to con- 
clude that the electron states describing the superfluid 
will no longer obey Bloch's theorem. 

Furthermore, we have recently proposed that the su- 
perconducting state should be understood as a (dis- 
torted) 'mirror image' of the normal state, where the 
'mirror' switches the role of the electron-ion and electron- 
electron interactions 3 . We argued that in the supercon- 
ducting state, the electron-electron interaction is priv- 
ileged over the electron-ion interaction and the many- 
body wavefunction optimizes the electron-electron in- 
teraction, with the electron-ion interaction being non- 
optimized. It is natural that such a state should be de- 
scribed by a single electron basis that does not satisfy 
Bloch's theorem. But what is the optimal basis to de- 
scribe such state? 

The new single-particle basis states that we propose 
here are defined by the relation 



'^k(r) = 



ipkjr) + ^tp-k(r) 
V2 



(5a) 



where i^Sk are the usual Bloch states for a given band. 
Note that {i^k} is a complete orthonormal set if and only 



2 



if {vk} is. The inverse relation is 



¥'k(r) 



V2 



(5b) 



Similarly wc define the electron creation operators in the 
new states cj^.^ in terms of the Bloch creation operators 

cL by 



Cka = 



^/2 



(6a) 



(6b) 



It is clear that our new basis functions do not satisfy 
Bloch's theorem, since 



(^k(r + R) 



e"'-"^k(r) + ze-^''-"-(/j_k(r) 



V2 



^ e''-''^k(r) 



(7) 



Nevertheless, I will argue that the new basis has very 
interesting properties and is relevant to the description 
of the superconducting state. 



II. NATURE OF THE WAVE FUNCTIONS 
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FIG. 1: Wave function amplitudes and charge densities for 
the standing waves associated with the operators c|^^ (full 
lines) and cl^ _^ (dashed lines) , assuming the Bloch states 
are plane waves. Note that the charge densities avoid each 
other reducing the Coulomb repulsion. 



Consider the case where the Bloch functions are well 
approximated by plane wave states 

^u(r) = (8) 

Then 

<Pk(r) = — ^cos(k • r - -) (9a) 

^-k(r) = ^cos{k -r + j). (9b) 

describing standing waves. The values of r for which 
|(^k(r)P is maximum satisfy 

TT 

k • r — - = niT (10a) 
4 

with n integer, and for those r values 

k.r+ J = (n+i)7r (10b) 

so that |(;5_k(r)P = 0. Consequently this will cause a 
strong reduction of the direct Coulomb repulsion between 
an electron in state <^k and another electron in state <^_k 
compared to the plane wave case where the density is uni- 
form, as illustrated in Fig. 1. Conversely, the in-phase 



modulation of the wavefunction amplitude for two elec- 
trons in the same state <^k will give a strong enhancement 
of the direct Coulomb repulsion in that case compared to 
the plane wave case. 

Similarly, consider a band described by the tight bind- 
ing approximation with a single orbital <pi{r) per site. 
The Bloch functions are 

Mr) = ^T.^''''''Mr) (11) 
and the new orbitals are given by 

^k(r) = ^ E ""^^ • - i)-^' W (12a) 

^ i 

(^_k(r) = i^Ecos(k.R, + J)</.,(r) (12b) 

It is clear that for small k the situation is the same as 
for the plane-wave case: the values of Ri for which the 
amplitude is maximum in Eq. (12a) coincide (or nearly 
coincide) with those where the amplitude in Eq. (12b) 
is zero. However, this changes when the wavevector k 
becomes an appreciable fraction of a reciprocal lattice 
vector G. We have 

TT 1 

cos(k • Ri =F 7) = ^(cos(k • R,) ± sm(k • Rj)) (13) 
4 \/2 



3 



so that the ampHtudes in Eq. (12) become the same when 
k — >■ G/2, since the second term in Eq. (13) vanishes. 
Thus, we expect the direct Coulomb repulsion between an 
electron in state i^k and one in state i^_k to be reduced, 
and that between two electrons in the same state i^k to 
be enhanced, when k is near the bottom of the band but 
not when k is near the top of the band. 



III. COULOMB REPULSION IN THE NEW 
BASIS 



The Coulomb repulsion between an electron in Bloch 
state k and one in Bloch state k' is 



<kt,k'4|Fe|kt,k'4>= 



U_ 
N 



(18) 



independent of k and k'. Instead, in the new basis 

< kt,k';|K|kt,k'i>= 

^(1 + sin{2k ■ Ri))(l + sin{2k.' ■ R;)) (19) 



The bare Coulomb interaction between opposite spin 
electrons in a plane wave basis is given by 

Vc = ^Y^ ^(q)4+qt4'-q4Ck'4Ckt (14a) 



kk'q 



^(q) 



(14b) 



with n the volume of the system. The Coulomb interac- 
tion between an electron of spin f in the plane wave state 
<^k and an electron of spin J, in ip'^ 



|kt,k' i>=4^4,jo> 

is simply 

< kt,kU|^c|kt,k'4>= 



(15a) 



/ dVdV|<^k(r)P^^|(pk'(r)P = ^V'(0)(15b) 



independent of k and k'. The Coulomb interaction be- 
tween an electron of spin t in state <^k and an electron 
of spin 4 in state 



kt,k'4>= 4t4'j0> 



(16a) 



is given by 



< kt,k'4|K|kt,k'4>=^no) 
+ ^l^(2fc)(-5k',-k + <5k',k) 



(16b) 



showing indeed a reduction for (k t, — k 4-) in the new 
basis as well as an increased repulsion for (k t, k J,), since 
the Coulomb matrix element V{2k) is always positive. 
Note also that the effect is largest for small k. 

Similarly consider the Coulomb interactions in the 
tight binding model. We will assume for simplicity that 
the only non-zero matrix element of the Coulomb inter- 
action in the tight binding model is the onsite Coulomb 
repulsion U: 

U = J dVdV|<^,(r)|2^^|<p,(r')|^ (17) 



showing in particular that the interaction is most reduced 
for k' = — k, and most enhanced for k' = k. However, for 
k and k' denoting states near the top of the band, ~ G/2, 
with G a reciprocal lattice vector, the effect disappears, 
in agreement with the discussion in the previous section. 

However, Eq. (19) appears to suggest that the 
Coulomb interactions are the same also for states around 
the half- filled band, e.g. k = G/4. This is an artifact 
caused by keeping only the on-site Coulomb matrix el- 
ement U. When considering other Coulomb matrix ele- 
ments also it is found that the Coulomb interaction be- 
tween states in the new basis and in the conventional 
basis become the same only when the states approach 
the top of the band and the wavevectors are connected 
by reciprocal lattice vectors. 

Next we consider the Coulomb interaction in the triplet 
and singlet (k t, — k J,) pairs. We have 



so that the triplet state in the new basis, 



t>-- 



'-kt^'-kj. 



'--kt'-kj. 



> 



(21) 



is identical to the triplet state in the conventional basis, 
\t >. The Coulomb interaction is 



< i\Vc\i>=< tlVclt >-- 



{V{Q) - V{2k)). (22) 



The singlet states are different however, since 



^kt'^-ki 



and we have 



'^-kt'^kj. 



kt'^-kj. 



< s\Vc\s >= 



< s\Vc\s >= -{V{0) + V{2k)) 



(23) 



(24a) 



(24b) 



so the Coulomb repulsion is smaller in the singlet state of 
the new basis than in the singlet state of the conventional 
basis. The average of Eqs. (22) and (24a) or (24b) yields 
Eq. (16b) or Eq. (15b) for k' = — k, since the state 
|k t, — k l> is a linear combination of singlet and triplet 
states with equal weight. 



4 



A similar result holds in the tight binding model, where 
we have (assuming only on-site interactions) 



< tWrlt >=< tWrlt >= 



(25a) 



21 J 

< 5|V;|S ^^cos22k-R, (25b) 



< sWcls > = 



2U 
"/V 



(25c) 



Again the average of Eq. (25a) and Eqs. (25b) or (25c) 
yields Eq. (19) or (18) for k' = -k. 



IV. RELEVANCE TO SUPERCONDUCTIVITY 

Our findings in the previous sections suggest that the 
new basis is relevant to superconductivity. The fact that 
the BCS wavefunction Eq. (4) singles out pairs (k f, — k | 
) could naturally be explained by the fact that the direct 
Coulomb repulsion between an f-spin electron in state 
(flk and a 4,-spin electron in (p-k is particularly small. 

Consider then as an alternative to the BCS wavefunc- 
tion a wavefunction of the form Eq. (4) but with the new 
states instead of the Bloch states: 



of the state ifk with the state (p-k being empty, and vice 
versa, and no amplitude for the states (^k and (^_k be- 
ing both singly occupied . Instead, the new state has 
amplitude for the states t^k and (p-k being both singly- 
occupied and no amplitude for the state ipk being doubly 
occupied with the state (^_k empty, nor vice versa. Sin- 
gle occupancy of both ip^ and ip-k minimizes the direct 
Coulomb repulsion between the electrons while double 
occupancy of either state maximizes the direct Coulomb 
repulsion. Consequently, we conclude that the new wave- 
function Eq. (26) will be strongly favored energetically 
over the BCS wavefunction Eq. (4) as far as the Coulomb 
interaction between electrons is concerned. 

In summary, in the normal state where states k and — k 
are both either occupied or empty, the new basis is com- 
pletely equivalent to the conventional one. However in a 
state of the BCS form that allows for partial occupation 
of pair states the new basis appears to be favorable. 

Furthermore, Eqs. (29), (16b) and (14b) indicate that 
the advantage of the new basis over the conventional one 
will be greatest if u^Vk 7^ in a region where k is small. 
Instead, for k near the edge of the Brillouin zone there 
is no advantage to the new basis according to Eq. (19). 
This suggests that in a superconducting state described 
by the wavefunction Eq. (26) the region where itk'i'k 7^ 
should occur near fc = 0, and the states with k near the 
Brillouin zone edge should be full (wk = l,Uk = 0) (to 
conserve the number of particles). This is precisely the 
scenario predicted in Ref. Q ■ 



Deep in the Fermi sea, Uk = U-k = and t^k = U-k = 1- 
Since 

(oUcUOicUAO = (4t^'-k,)(clkt4,) (27) 

deep in the Fermi sea the new wavefunction is the same 
as the BCS wavefunction. This also means that in the 
normal state (wk^k = for all k) the wavefunctions Eq. 
(4) and Eq. (26) are identical. However, in the region 
where Mk^k 7^ there is a difference. We have for (k f 
, — k I) and (— k f, k \.) pairs in the BCS wavefunction 

l*>BCS = ('^k + Wkc[tC-ki)("-k + ^^-kC^_kt4i) 



^k (Ckt<^-k4,) ('^-kt^kj.) 



- *"ki'k((4T4;) + (Siktclki)) (28) 

(assuming as usual Uk — u_k, Wk — ''^-k)- Instead, in the 
new wavefunction Eq. (26) 

\'^>new^ = (wk + Wk'5ktctki)(^^-k + ^'-k6lkt4i) 
= '4 + ^'k(4tC-ki)(clkT4i) 
+ UkVk{{£l^£l^^) + (sLkt4;)) (29) 

Eqs. (28) and (29) differ in the last line. The BCS wave- 
function has an amplitude for having double occupancy 



V. 4-ELECTRON STATE 

We consider in the following the plane wave case only. 
Let us consider further the two candidate wavefunctions 
Eqs. (4) and (26) for the case of only two pairs, in states 
(k t, -k i) and (-k t, k |). The Coulomb energy in the 
BCS case is 



4^(0), ,4 ^ 2(1/(0) + y(2fc)), ,2, ,2 
< Vc >BCS= —^\vk\ + ' kkl bkl 



n 



and with the new wavefunction Eq. (26) it is 



4^(0), ,4 



2^0), ,2, ,2 



(30a) 



(30b) 



Since the Coulomb repulsion V{2k) is positive for all fc, 
the new state has clearly lower Coulomb energy than the 
BCS state. 

However in adding the second electron pair we have 
lost the lowering of energy proportional to V{2k) that 
occurred for the single pair (Eq. (16b)) due to interfer- 
ence effects. How can it be restored? In deriving Eq. 
(30), we assumed Wk = W-k, ^k = I'-k- Not making that 
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assumption we obtain instead 

< V,>Bcs=^i\v^.\' + \v-^.\Y 



V(2fc) 

n 



Wk^^kU-ki'* k + W-k^-kWki'k) (31a) 



and with the new wavefunction Eq. (26) 

+ ^(i^kpbkp+Kknt'-kn 

V{2k) 



2n 



The last term in Eq. (31b) lowers the energy if 
Uk u-k 



(31b) 



(32) 



that is, if the occupation of (k — k l) is different from 
that of (— k t, k J,) . Thus, it favors a spontaneous break- 
ing of parity and the existence of a spin current [9] . 

In earlier work we have presented several arguments 
in favor of the suggestion that superconductors possess 
a spontaneous spin current in their ground state 0, 0- 
0|. The direction of the spin current is determined by 
the spin-orbit interaction in the presence of an outward- 
pointing electric field predicted by the theory of hole 
superconductivity [Tl'l together with the sample geome- 
try. Namely, near a surface with outward normal fi the 
member of the Cooper pair of spin a has larger velocity 
in direction n x . Also, C.R. Hu has argued quite gen- 
erally that the structure of the BCS state is compatible 
with the existence of a spin current [l2| . 



VI. COULOMB INTERACTION IN THE PAIR 
WAVEFUNCTIONS 

The average of the Coulomb interaction in the BCS 
wavefunction Eq. (4) is given by 



< Vc >BCS = 



kk' k 

(33) 

and in the new wavefunction Eq. (26) by 

<V,>nen, = ^ ^ < W<JV^(k' - k) - ^ (k' + k)] 
kk' 

k 

V^(0) V~~^, ^ * I |2| \2\\ 

+ Z^("kWkU_kW_k + |Wk| \Vk\ )) 

k 

+ ^E^(2fc)KP[Kkp-|«kP) (34) 



For a clearer comparison we separate self-energy terms 
in the double sums and obtain for the BCS case 

< Vc >BCS^ E UkWkUk'Wk'^(k- k') 



^(0) 

n 



(E 



k^k' 



El^kH^kP (35) 



and for the new wavefunction 



^ E KvkUk<>[Vik' -k) - V{k' + k)] 



k^^k' 
k^-k' 



^(0) 

n 



^ ^ V{2k)\ukV-u - u-kWkl 



(36) 



The BCS expression Eq. (35) has the same form 
whether or not parity is broken. Instead, in the ex- 
pression Eq. (36) both the first and last term vanish 
identically if parity is unbroken, and the average of the 
Coulomb interaction (as well as of any other interaction) 
gives a self-energy which is identical to that of the BCS 
case (2nd lines in Eqs. (35) and (36)). 

Thus we conclude that there is a stark difference be- 
tween the BCS wavefunction and the new wavefunction. 
The BCS wavefunction will give rise to superconductiv- 
ity if the first term in Eq. (35) is negative on average's]. 
This will be the case if the interaction V{k — k') in addi- 
tion to the repulsive Coulomb interaction contains neg- 
ative (attractive) terms arising e.g. from the electron- 
phonon interaction, and if the attractive terms are suf- 
ficiently large that they can overwhelm the repulsive 
Coulomb terms. Of course the frequency dependence 
of the electron-phonon interaction provides additional 
help by reducing the effective strength of the Coulomb 
interaction [l3|. 

Instead, with the new wavefunction an attractive inter- 
action (whether instantaneous or retarded) will not help 
if parity is unbroken, since it will cancel identically in the 
first term in Eq. (36). The direct Coulomb repulsion is 
however also eliminated in the first term in Eq. (36) in 
contrast to the first term in Eq. (35) if parity is unbro- 
ken. Furthermore, if parity is broken, the last term in 
Eq. (36) provides a lowering of energy in the supercon- 
ducting state for a repulsive interaction^ and particularly 
if the region where Uk^k 7^ occurs near fc = 0. 



VII. DISCUSSION 

It is natural to try to describe the superconducting 
state using electronic states that do not optimize the 
electron-ion interaction and many such attempts have 



6 



been tried in the past. For example, Froehlich^Ti'] pro- 
posed a state where a sheU of occupied states at the top of 
the normal Fermi distribution is moved outward, leaving 
behind an empty shell. This has a certain resemblance to 
the scenario proposed in Ref. @. Even in conventional 
BCS theory, the normal state Fermi surface is 'blurred' 
and Bloch states that are empty in the normal state (i.e. 
unfavorable for the electron-ion interaction) become par- 
tially occupied in the BCS state. However we are not 
aware of earlier attempts to describe the superconduct- 
ing state using single particle states that violate Bloch's 
theorem as we pro pose in this paper, except for the early 
work of Kronig[l5i and Heisenberg|la| that proposed the- 
ories of superconductivity based on an electronic 'lattice' 
that had no relation to the underlying ionic lattice. In 
addition, Ref. [17] used sine functions to take into account 
the effect of sufaces and impurities. 

As discussed here, the basis introduced in this paper 
is consistent with the proposal of ref. ^ that in the su- 
perconducting state the states near the top of the band / 
edge of the Brillouin zone become full and those near the 
bottom of the band / center of the zone, that resemble 
free-electron-like states, become empty or partially occu- 
pied (Fig. 6 of ref. |[6j]). This is because the new basis 
gives rise to energy lowering compared with the conven- 
tional basis for small k (near the bottom of the band). 



while it has no advantage over the conventional basis for 
k near the edge of the Brillouin zone (states near the 
top of the band). This in turn is consistent with the 
observation that properties of the superconducting state 
are 'universal', because they arise from long wavelength 
states near the bottom of the band that are insensitive to 
details of the electron-ion potential. In order for a sys- 
tem to be able to push the electronic occupation to the 
top of the band (Fig. 6 of ref. @) without undue energy 
cost requires the band to be almost full in the first place 
(hole conduction in the normal state), and is favored by 
a low curvature of the band near the top (large effective 
mass), a narrow band and a small ionic charge Z, which 
are all requirements that we discussed earlier in connec- 
tion with the theory of hole superconductivity but based 
on different arguments. 

We also found here that with the new basis the energy 
is lowered if parity is broken and a spin current exists, 
consistent with our earlier predictions [H, |^. In a small 
cluster, a rigid electronic state where parity is broken 
should lead to an electric dipole moment as observed in 
the remarkable experiments of de Heer and coworkers [Tsf. 

The new basis proposed here may be of interest also for 
other unconventional theories of superconductivity be- 
sides the theory of hole superconductivity. 



[1] F. Bloch, Z. Physik 52, 555-600 (1929). 
[2] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 - A1138 
(1965). 

[3] J. Bardeen, L. N. Cooper , and J. R. Schrieffer, Phys. 

Rev. 108, 1175(1957). 
[4] J.E. Hirsch, Phys.Rev. B71, 104522 (2005). 
[5] J.E. Hirsch, Europhys. Lett. 81, 67003 (2008); J. Phys. 

Cond. Matt. 20, 235233 (2008). 
[6] J.E. Hirsch, Int. J. Mod. Phys. 23, 3035 (2009). 
[7] J.E. Hirsch, Phys.Lett.A 309, 457 (2003). 
[8] J.E. Hirsch, Phys.Lett.A 315, 474 (2003). 
[9] J.E. Hirsch, Phys.Rev. B71, 184521 (2005). 
[10] J.E. Hirsch and F. Marsiglio, Phys.Rev. B39, 11515 

(1989). 



[11] J.E. Hirsch, Phys.Rev. B68, 184502 (2003). 
[12] C.R. Hu, arXiv: 0711.3020 (2007). 

[13] P. Morel and P.W. Anderson, Phys. Rev. 125, 1263 
(1962). 

[14] H. FroeUch, Phys. Rev. 79, 845 (1950). 
[15] R. de L. Kronig, Z. Phys. 78, 744 (1932). 
[16] W. Heisenberg, Zeitschrift fur Naturforschung 2a, 185 
(1947). 

[17] K. Tanaka and F. Marsiglio, Phys.Rev. B 62, 5345 
(2000). 

[18] S. Yin, S. Xu, A. Liang, J. Bowlan, R. More and W.A. de 
Heer, Jour, of Superconductivity and Novel Magnetism 
21, 265 (2008) and references therein. 



